Bragg spectroscopy with an accelerating Bose-Einstein condensate 
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We present the results of Bragg spectroscopy performed on an accelerating Bose-Einstein conden- 
sate. The Bose condensate undergoes circular micro-motion in a magnetic TOP trap and the effect 
of this motion on the Bragg spectrum is analyzed. A simple frequency modulation model is used 
to interpret the observed complex structure, and broadening effects are considered using numerical 
solutions to the Gross-Pitaevskii equation. 
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PACS numbers: 03.75.Kk, 03.75.Nt, 32.80, 32.80.Cy 



I. INTRODUCTION 

The experimental realization of a dilute Bose-Einstein 
condensate (BEC), a quantum-degenerate atomic Bose 
gas, has led to a wide range of experiments in atom optics 
(see, for example, Ref. [3,I3)- The macroscopically oc- 
cupied ground state exhibits large-scale coherence, which 
can be exploited in atom optical techniques. One such 
technique is Bragg scattering, which has been applied 
to make precise measurements of condensate momentum 
0, 01 , investigate excitations , measure the phase co- 
herence length [^jexcite phonons and to implement 
four-wave mixing [3| and a Mach-Zehnder interferometer 
0. More recently, Blakie et al. showed that the veloc- 
ity sensitivity of Bragg scattering can provide a means 
of detecting superfluid flow associated with condensates 
in a vortex state In all this work, the initial Bose 

condensates have been well described as wave-functions 
with a stationary center of mass (COM), making them 
an ideal source for atom-optic experiments using Bragg 
scattering. 

In conventional Bragg spectroscopy the spectrum is 
given by the spectral response function of the conden- 
sate which is related to the momentum distribution 
and under certain conditions is a good measure of this 
distribution. In this work, we show that for an ac- 
celerating condensate, the spectrum is no longer a di- 
rect representation of the momentum distribution. We 
investigate Bragg spectroscopy for the case of a Bose 
condensate undergoing circular micro-motion in a Time- 
averaged Orbiting Potential (TOP) trap [H El. Al- 
though the condensate micro-motion is small in ampli- 
tude, we demonstrate that the COM acceleration can be 
large enough to dramatically alter the Bragg spectrum. 
Our experimental results are compared to a theoretical 
spectrum for the corresponding stationary condensate, 
illustrating that the effect of the acceleration is to signif- 
icantly broaden the Bragg spectrum and introduce com- 
plex structure. We describe two simple models to explain 
the underlying physics of the observed behavior, and 
present numerical simulations using the one-dimensional 



Gross-Pitaevskii equation (CPE) which confirm our in- 
terpretation. 



II. CONDENSATE MICRO-MOTION AND 
BRAGG SPECTROSCOPY 

In this section we summarize aspects of micro-motion 
in a TOP trap that are relevant to its affect on Bragg 
spectroscopy, and present two simple models for describ- 
ing the modified Bragg scattering process. 

The instantaneous magnetic field for our TOP trap is, 
as given in Ref. p^. 



B = {B'^x+Bt cos{ujTt))i+{B'^y' 

(1) 

where B'^ is the quadrupole field gradient along the ra- 
dial direction, Bt is the magnitude of the rotating bias 
field, and ujt is the angular frequency of the the bias 
field rotation. Petrich et al. showed that, for time 
scales longer than the rotation period, the atoms experi- 
ence only the leading terms in the magnitude of the time 
averaged magnetic field so that 
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where r is the radial coordinate. For many experiments, 
this time-average is an appropriate description, but for 
Bragg spectroscopy this is not always the case [T^- The 
time-varying force on the condensate results in COM os- 
cillation or micro-motion (as pointed out by Petrich et 
al. 101), to which Bragg spectroscopy is sensitive. The 
origin of this micro-motion can be seen by considering 
the instantaneous potential of the TOP trap 



U{x,y,z,t) = ^|B| 

= MS;[r2-}-a;2-hy2+4^2 

+2rQ{x coii{ujTt) + y s\n{ujTt))]^ , (3) 

where ro — Bt / B'^ is the so-called "circle of death" 
radius (the radius where the B = Q point rotates in 
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the z = plane) and the time-varying force is given by 
F{x,y,z,t) = —VU. Typically, in the region of the con- 
densate in the trap, tq ^ x,y and we can approximate 
the magnitude of the force along the x-axis (for example) 
as 

F{t)^-fiB'^cos{iJTt). (4) 

Ignoring condensate mean-field effects (since the trap- 
ping force is large) we find that the classical solution for 
the magnitude of the time- varying COM momentum and 
position are given by 

p{t) = -—^sm{uJTt) (5) 

OJt 

x{t) — 5-cos(wTi), (6) 

where m is the atomic mass and we have assumed for sim- 
plicity p(0) — 0, and a;(0) = [iB'^jmio^. The condensate 
COM along the a;-axis performs simple harmonic motion, 
and from trap symmetry we deduce that the condensate 
performs circular motion in the plane of the rotating bias 
field. This result was first obtained and experimentally 
confirmed by Miiller et al. 11] . 

To discuss the implications of centripetal accelera- 
tion for Bragg spectroscopy we first briefly summa- 
rize Bragg scattering 14] . A moving optical grat- 
ing is applied to the condensate and atoms can ab- 
sorb n photons from one beam (wi , ki ) and emit 
n photons, via stimulated emission, into the other 
(uj^i'k.-i) where ]ki] = ]k2l = k = 27r/A, and the detuning 
is S = coi — UJ2- This results in the transfer of momentum 
h\ci\ = nPrecoii = 2nhk siii{9 / 2) and energy hnd, where 
q = ki — k2 , and 9 is the angle between the two beams 
forming the grating. This process will be resonant for 
non-interacting stationary atoms if the final kinetic en- 
ergy of the atom is equal to the energy difference between 
the Bragg beams i.e. Un5 = {nPrecoiif /2m. 

When calculating the resonance condition for atoms 
which are initially moving in the lab frame an additional 
term arises, which can be attributed to the Doppler ef- 
fect, and which we will call the Doppler term TiP ■ q/m 
(where P is the condensate momentum). In the micro- 
motion case we have P = p{t)i giving a Doppler term 
which results in a new time-dependent Bragg detuning 
condition. For first order Bragg scattering (n = 1) this 
condition is 



5 < Precoii/'^''^^ resonance only occurs if sin{uJTt) > 0. 
Therefore, if the Bragg pulse length, Tg, is not an inte- 
ger multiple of the TOP period the Bragg spectrum will 
be asymmetric. We avoid this problem by having Tb 
equal to eight TOP rotations. 

Alternatively, we can picture the Bragg scattering pro- 
cess in the rest frame of the condensate and think of 
the micro-motion as modifying the optical potential. In 
the absence of micro-motion the effective optical poten- 
tial from two counter-propagating, pulsed, laser beams is 
given by 

Voptix,t) = hV{t)cos{\q\x-6t), (8) 

where the amplitude hV{t) is twice the AC stark shift 
at the grating maximum. In this case, the potential 
has a spectrum centered at the frequency difference be- 
tween the two Bragg beams, with an rms half-width 

Ai^Tu ^ i/(2rs). 

The effect of the micro-motion can be incorporated into 
the optical potential by transforming to the stationary 
frame of the condensate. The resulting time-dependent 
frequency shift of the Bragg resonance condition in the 
laboratory frame produces a frequency modulated (FM) 
optical potential in the condensates rest frame, as pointed 
out by Cristiani et al. The modified potential be- 

comes 

Voptix^t) = hV{t) cos{\q\x -St - \q\Acos{ujTt)), (9) 

where the bar denotes the COM frame of the condensate 
and A = fiB'^/muj^. This frequency modulation gener- 
ates side-bands in the frequency spectrum of the optical 
potential, yielding multiple first-order Bragg resonances. 

Bragg spectroscopy involves scanning the frequency 
difference 6 through resonance and measuring the out- 
coupled fraction. To resolve features in the conden- 
sate's Bragg spectrum requires the frequency width of 
the Bragg light to be narrower than the spectral width 
of the feature. Mechanisms that contribute to the con- 
densate momentum width 0, Q are Doppler broaden- 
ing Ai'D = y/ {21/8){Precoii/2TrmR), where R is the con- 
densate radius, and mean-field (collisional) broadening 
Az/„ = \/{8/M7){njnax/h), where fi^ax is the chemical 
potential at the peak density. In addition to these ef- 
fects and the spectral broadening associated with finite 
pulse time, there is also a power broadening contribution 
Afp ~ kV/T:\q\ J^]. 
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(7) 



Note that when S is greater than Precou/'^^^ reso- 
nance condition can only be satisfied for su\{u>j't) < 0, 
so that Bragg scattering is resonant only for some 
instant during half of the TOP cycle. Similarly, if 



III. EXPERIMENTAL DETAILS 

Our apparatus for producing Bose condensates is de- 
scribed in ^10], but with some minor modifications. 
Briefiy, we use a double magneto-optical trap (MOT) 
to cool and trap *^Rb atoms, which is driven by 
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injection-seeded diode lasers and with continuous load- 
ing of the low-pressure MOT. The laser cooled sample 
is transferred to a magnetic TOP trap and evaporatively 
cooled to produce a condensate containing approximately 
2 X 10^ atoms in the \F — 2,mp — 2) hyperfine ground 
state. 

Following condensate formation, the magnetic trap was 
adiabatically relaxed (by reducing the quadrupole field 
gradient and increasing the bias field over 200 ms) to 
trapping frequencies of uJradiai = i^z/n/S = 27r x 18 Hz. 
For our bias field rotation frequency lot — 2hx 2.78 kHz, 
the corresponding tangential micro-motion velocity is 
3.3 mm s~^ (compared to the two-photon recoil veloc- 
ity of 11.8 mm s~^). Although the amplitude of the mo- 
tion is only 190 nm, the high relative velocity means that 
the effect of the micro-motion on the Bragg spectrum is 
substantial. 

Counter-propagating Bragg beams were generated us- 
ing a Coherent MBR Ti:sapphire laser, the output of 
which was split and passed through two acousto-optic 
modulators (AOMs) to provide individual frequency con- 
trol. We used a detuning A ~ 27r x 4.49 GHz above the 
F = 2 to F' ~ 3 transition, giving a sponta- 

neous scattering rate of approximately 5 (i.e. neg- 
ligible for our purposes). Both beams were aligned in 
the plane of the rotating bias field. A pulse length of 
Tb — 2.88 ms was chosen for a spectral width Ai/Tb com- 
parable to the width of an equivalent stationary con- 
densate with our parameters. The pulse length is also 
shorter than a quarter trap period so that the Bragg 
scattered fraction becomes well separated from the un- 
scattered component in time-of- flight. The intensity was 
set between 0.26(5) to 0.6(1) mW cm~^. To minimize 
the effect of timing differences between consecutive mea- 
surements we synchronized the Bragg pulse with the ro- 
tation of the bias field. At the end of the pulse, the 
condensate was released from the trap and imaged by 
resonant absorption after 9 ms of free expansion. Bragg 
spectroscopy was performed by measuring the first-order 
scattered fraction as a function of the frequency differ- 
ence between the two Bragg beams. 



IV. RESULTS AND DISCUSSION 

For a stationary Bose condensate confined in a static 
harmonic potential with our experimental parameters, 
the Bragg spectrum has a simple form (essentially Gaus- 
sian) with a narrow width. However, the presence 
of micro-motion introduces substantial broadening and 
complex structure. In Fig. ^ we show the Bragg scat- 
tered data obtained for a condensate undergoing cen- 
tripetal acceleration in the plane of the Bragg beams. 
Also shown is a theoretical spectrum for our condensate 
parameters, ignoring mean-field effects (which are small 
in this case) and assuming no micro-motion. The rms 
half- width of the theoretical curve is 0.56 kHz, which is 
consistent with our measurements made in time-of-flight 
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FIG. 1: Bragg scattered data for a Bose condensate accel- 
erating due to micro-motion in a TOP trap (crosses), to- 
gether with a numerical simulation using the linear ID Gross- 
Pitaevskii equation assuming no micro-motion (solid line). 
Results are for = 90 G cm~^ (corresponding to a max- 
imum COM velocity of 3.3 mm s~^) and an intensity of 
0.26(5) mW cm~^ (corresponding to a power broadening con- 
tribution of 0.22(4) kHz). 



where micro-motion is absent. 

To explain the complex structure in the experimental 
spectrum, we now present results of the simple FM pic- 
ture described in the previous section. We obtain the 
frequency components present in the Bragg light by per- 
forming a Fourier transform on the modified optical po- 
tential given by Eq. The result is a spectrum for the 
Bragg light which has a carrier at the frequency difference 
between the Bragg beams (6), and side-bands (spaced by 
lut) with relative amplitudes determined by Bessel func- 
tions of the FM spectrum ^^l- In this case the modula- 
tion index, upon which the Bessel functions are strictly 
dependent, is just the maximum Doppler shift divided 
by the modulation frequency and can therefore be ad- 
justed by simply changing the quadrupole field gradient. 
In Fig. 13 the experimental data from Fig. ^ is overlaid 
with the corresponding frequency spectrum of the mod- 
ified optical potential Vopt (scaled by an arbitrary fac- 
tor). Since there are a number of side-bands with sig- 
nificant amplitude, when we sweep the carrier frequency 
(S) to generate the Bragg spectrum we bring successive 
side-bands into resonance with the condensate. The ex- 
act form of the Bragg spectrum is then a combination 
of the complex spectral character of the light and the 
spectral response of the stationary condensate. When 
the Bragg spectrum of the stationary condensate is nar- 
row compared to the bias field rotation frequency (as in 
our case), the Bragg spectrum with micro-motion has 
peaks at frequencies corresponding to the locations of 
the side-bands. To within the spectral uncertainty of our 
Bragg experiments, approximately ±1 kHz, this behavior 
is shown in the data of Fig. [3 Although we cannot re- 
solve the predicted individual spectral peaks, the overall 
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FIG. 2: The Bragg scattered data from Fig. Q (crosses) , to- 
gether with the spectrum of the Bragg hght (dotted line) and 
the corresponding numerical simulation (solid line) using the 
linear ID Gross-Pitaevskii equation. The relative heights of 
the peaks in the FM spectrum are determined by the modu- 
lation index and the width of each peak is determined by the 
finite pulse length Tb- Each data point has an uncertainty 
in the scattered fraction of ± 0.1. The dashed line indicates 
the minimum detectable fraction for this data and below this 
limit we plot the data as having zero scattered fraction. 



range of frequencies for which we get Bragg scattering 
matches the range over which there is significant side- 
band intensity. Our spectral resolution is limited by a 
number of physical processes, including the mechanical 
stability of the optics and the electrical stability of the 
currents driving the magnetic trap. From a heterodyne 
measurement using the two Bragg beams we determined 
that the uncertainty due to the mechanical stability of 
the optics (which Doppler shifts the light) is approxi- 
mately ±0.4 kHz, whereas the contribution to the uncer- 
tainty from current stability is approximately ±0.6 kHz. 
The main effect of current instability is to put noise on 
the quadrupole gradient, which in turn puts noise on 
the time-varying COM momentum, as can be seen from 
Eq. (0). The uncertainties in our measurements of the 
fraction of scattered atoms were determined by making 
repetitive measurements with the same parameters, and 
can vary across the spectrum. The significant contri- 
butions arise from shot-to-shot condensate variation and 
the minimum scattered fraction we can measure with our 
imaging system. 

At low intensity the simple FM picture predicts some 
important features such as the frequencies and relative 
heights of the peaks, but the widths of each peak only 
include a contribution from the finite pulse time. To 
obtain a more accurate prediction of the experimental 
results we include the dynamic term of the TOP trap in 
a collisionless one-dimensional Gross-Pitaevskii equation 
simulation, using a one-dimensional version of Eq. ^ as 
our magnetic potential. Results of this simulation are 
shown as the solid line in Fig. [3 confirming aspects of 
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FIG. 3: Power-broadened Bragg scattered data (crosses) and 
the corresponding ID linear GPE simulation (solid line). Ex- 
perimental parameters are = 90 G cm~^ and a intensity of 
0.6(1) mW cm^^ (corresponding to a power broadened contri- 
bution of 0.5(1) kHz). Each data point has an uncertainty in 
the scattered fraction of between ± 0.05 to ± 0.09 depending 
on 5/2n. 



the FM picture and providing better agreement with the 
data. 

As an additional test of our interpretation wc now 
consider a power broadened case (i.e. a higher value for 
V{t)). This alleviates the need to resolve very narrow 
spectral features and gives larger (more easily measured) 
Bragg scattered fractions near the central region of the 
spectrum. Fig. |31 shows a power broadened Bragg spec- 
trum together with the corresponding numerical simu- 
lation. The contribution to the condensate width from 
power broadening is 0.5(1) kHz, compared to 0.26 kHz 
from Doppler broadening (the dominant condensate con- 
tribution). Within the limits of our experimental uncer- 
tainty, the agreement between the data and the simula- 
tion is reasonable. The overall width of the spectrum 
is well predicted, and the variation in the Bragg scat- 
tered fraction is reduced, as expected. The small differ- 
ence between the center frequencies of the GPE simula- 
tions and experimental data spectra is likely to be due to 
our linear one-dimensional modeling of a nonlinear three- 
dimensional experiment 0, ITsl IT^ . The slight asymme- 
try in the numerical spectrum (which we cannot resolve 
in the experimental data) becomes apparent when trans- 
forming to the condensate rest frame through a damped 
drift term . Additional asymmetry can occur if the 
Bragg pulse length does not match an integer multiple 
of the TOP trap rotation period, but this is not the case 
for the results presented here. 

A process not included in our linear GPE simulations, 
but which can significantly alter the Bragg spectrum , 
is the mean field coUisional interaction in the condensate. 
A detailed theoretical treatment is currently being pre- 
pared [T^ . 
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V. CONCLUSIONS 

We have presented experimental results and a simple 
theoretical analysis of Bragg spectroscopy performed on 
an accelerating Bose-Einstein condensate. This work ex- 
tends Bragg scattering experiments beyond those involv- 
ing Bose condensates with a stationary (or constant ve- 
locity) center of mass. Using condensate micro-motion 
in a magnetic TOP trap, we have shown that for circular 
center of mass motion the Bragg spectrum is modified sig- 
nificantly from a simple momentum interpretation. An 
analysis which treats the angular acceleration as equiva- 
lent to frequency modulating the optical grating provides 
a useful interpretation for the complex structure we ob- 
serve in the Bragg spectra. Niimcrical simulations us- 
ing the one-dimensional linear Gross-Pitaevskii equation 
extend the analysis to include other broadening mech- 
anisms. Although we are not able to resolve the nar- 



row features predicted by our models, the overall spectral 
width and scattered fractions are in quantitative agree- 
ment with the theoretical results, confirming our simple 
physical interpretation. However, a more comprehensive 
theoretical description of Bragg spectroscopy with accel- 
erating Bose condensates is needed to characterize the 
efi'ects of collisions and collective excitations in higher 
dimensions. 
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